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Abstract 

The detection of a "Cold Spot" in the CMB sky could be explained by the presence of an anoma- 
lously large spherical underdense region (with radius of a few hundreds Mpc//i) located between us 
and the Last Scattering Surface. Modeling such an underdensity with an LTB metric, we investigate 
whether it could produce significant signals on the CMB power spectrum and bispectrum, via the 
Rees-Sciama effect. We find that this leads to a bump on the power spectrum, that corresponds to 
an 0(5% — 25%) correction at multipoles 5 < I < 50; in the cosmological fits, this would modify 
the x 2 by an amount of order unity. We also find that the signal should be visible in the bispec- 
trum coefficients with a signal-to-noise S/N ~ 0(1 — 10), localized at 10 < £ < 40. Such a signal 
would lead to an overestimation of the primordial /nl by an amount A/jvz, — 1 for WMAP and by 
Af NL ~ 0.1 for Planck. 

PACS numbers: 98.80.Cq, 98.80.Es, 98.65.Dx, 90.70.Vc 



I. INTRODUCTION 

The recent WMAP [l[ experiment has measured with great accuracy the anisotropies of the Cosmic 
Microwave Background (CMB), whose features are in good agreement with the expectations from inflation 
of a Gaussian spectrum of adiabatic fluctuations, fully described by a nearly scale- invariant power spectrum. 
However, it has been pointed out by several authors that the data contain various unexpected features. 
Some of them are localized at very large angular scales, such as the low quadrupolc, the alignment of the 
low multipoles (for a review see Q and references therein) and the power asymmetry between the northern 
and southern hemispheres [|| . Another anomaly is the presence of the so-called Cold Spot [HSH, 

which 

is a large circular region on an angular scale of about 10° that appears to be anomalously cold: the 
probability that such a pattern would appear from Gaussian primordial fluctuations is estimated to be 
about 1.8% d, [H, ||. So, while this could still be a statistical fluke, some authors have put forward the 
idea that it could instead be due to an anomalously large spherical underdense region of some unknown 
origin, on the line of sight between us and the Last Scattering Surface (LSS) @, S|- We may also remind 
that [lj| has claimed that, looking at the direction of the Cold Spot in the Extragalactic Radio Sources 
(the NVSS survey), an underdense region is visible at redshift z ~ 1 (see, however [ll[ for a paper that 
challenges this claim). Another motivation for studying such objects is that an underdense region of two 
or three hundreds of Mpc/ h could be enough to give an acceptable fit to the Supernova data and the CMB 
without Dark Energy [12j, if we happen to live near its centre. 

In this paper we also take the point of view that the Cold Spot could be due to such an underdense 
region, customarily denoted as a "Void" . By modeling it through an inhomogencous Lemaitrc-Tolman- 
Bondi (LTB) metric (which also requires an overdense compensating shell), we compute the impact of such a 
Void on some observational quantities, focusing on the statistical properties of the CMB: in particular on the 
two-point correlation function (power spectrum) and the three-point correlation function (the bispectrum). 
This is interesting for the following reasons. First, if there is such a Void, does the power spectrum get a 
sizable correction? And, if yes, to what extent the estimation of the cosmological parameters is affected? 
Second, if there is such a Void, could it be detectable also in the bispectrum? Third, does the presence of 
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such an object interfere with the measurement of a primordial non-gaussianity, which constitutes a very 
important piece of information in order to distinguish between models of inflation? 

It is well known that, passing through a Void, photons suffer some blue-shift due to the fact that the 
gravitational potential is not exactly constant in time - the so-called Rees-Sciama (RS) effect By RS 
effect we mean the one associated to the variation of the gravitational potential from non-linear effects. 
There can be an effect already at the linear level, usually referred to as the Integrated Sachs- Wolfe (ISW) 
effect, which however would vanish in a matter dominated flat Universe. The ISW effect would be significant 
only when a Dark Energy component becomes dominant with respect to matter. Hence, here we focus our 
attention only on the RS effect, which is always present, and briefly comment on the extension to a ACDM 
Universe later on. Note that there is a second physical effect due to a Void on the line of sight, in addition 
to the RS blue-shift of the photons: the lensing of the primordial perturbations, which will be analyzed in 
detail in a companion paper [lil ]. 

It is also known that the RS effect scales as the third power of the comoving radius of the Void L 
times the present value of the Hubble parameter Hq, AT/T oc —(LHq) 3 . For Void sizes compatible with 
the expectations from the usual structure formation scenarios, the RS effect happens to be suppressed 
with respect to the primordial temperature fluctuations (l5j . In order to produce a signal comparable 
to the measured temperature anisotropy of \AT/T\ ~ 10~ 5 , the Void should be of very large size, i.e. 
L ~ (200 — 300)Mpc//i: this would be at odds with the standard scenario of structure formation from 
Gaussian primordial fluctuations (at more than 10cr @]). We may also mention that according to 0, E3] 
there are many localized regions in the sky (both underdense and overdense) of about 100Mpc//t, which 
would be responsible for the correlations between the CMB and the Large Scale Structures: as [l8[ has 
recently stressed the existence of these regions is already at odds with the usual structure formation scenario. 
Even though we do not address in this paper the issue of the primordial origin of such large objects, we 
mention some possibilities. For example, one could consider certain models of inflation, such as the ones of 
the "extended" type [H, HO, HH , with the possibility of tunneling events and nuclcation of bubbles, which 
would appear as Voids in the sky today. Alternatively, one could also imagine the presence of non-Gaussian 
features in the primordial fluctuations which would seed a large Void or even non-conventional structure 
formation histories in the late-time Universe which could enhance the probability of having large Voids. 

In the case of an explanation of the Voids via tunneling events, it is interesting to note that if a nucleation 
process has small probability (per unit volume and time) compared to the Hubble rate (at some time during 
inflation) then the number of anomalous Voids could well be very small, such as having just one or a few 
of them in the present observable Universe. In addition, since the LSS is a rather thin shell whose volume 
is much smaller than the total volume inside it, if there are only few Voids it is more likely that they are 
located along the line of sight, rather than at the LSS 1 . This is an important point because a Void at the 
LSS would have a huge impact on the CMB and, as a consequence, its size would be strongly constrained: 
L < 4Mpc//i [13, [Hj]. On the contrary, Voids on the line of sight are not subject to the latter constraint. 

Note also that other authors have put forward alternative ideas, such as the idea that the Cold Spot 
could be due to a topological defect, in particular a cosmic texture [23l |. As we are going to discuss, some 
of our considerations apply in a similar fashion if the origin of the Cold Spot is a texture, rather than a 
Void. 

The paper is organized as follows. In sect. |H] wc calculate the RS-induccd temperature profile of the 
CMB photons passing through a large Void with the physical characteristics of the Cold Spot. The impact 
on the power spectrum is discussed in sect. IIII[ while sect. IIVI deals with the impact on the bispectrum 
and on the overestimation of /jvx that would be done by neglecting the presence of such a Void. Finally, 
in sect. E]we extend some of the considerations to the case of several objects in the sky. In fact, it is 
conceivable that if there is one such structure, there may be other ones, maybe of smaller size and less 
visible - see e.g. [Till. [l7j for some other candidates in the CMB sky. Our conclusions are drawn in scct. lVIl 



It is possible, also, to imagine that some primordial process could generate a coherent spherical region on the LSS surface, 
with small density contrast, on top of the primordial Gaussian spectrum, and which would produce a 10~ 5 fluctuation. 
However, we do not consider here such possibilities. 
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II. A VOID IN THE LINE OF SIGHT 



As anticipated in the introduction, we would like to consider the following cosmological configuration: 
an observer looking at the CMB through one spherical Void located at comoving distance D from him. We 
assume that the object does not intersect the LSS and that the observer is not inside it. The observer re- 
ceives from the LSS the CMB photons, whose fluctuations we assume to be adiabatic, nearly scale-invariant 
and Gaussian (as those generated e.g. by the usual inflationary mechanism). Given this configuration, 
the observer detects one particular realization of the primordial inflation-generated perturbations on the 
LSS plus the secondary effects due to this anomalous structure located in the line of sight. Our aim is to 
give the theoretical prediction for the two-point and three-point correlation functions, in order to compare 
them with the observations. 

As a first step we consider that the Universe and the LSS are just completely homogeneous, and study 
what kind of profile we should see in the sky, because of the existence of one inhomogeneous Void located 
between us and the LSS: we denote the temperature fluctuation obtained in this way as AT^ RS ^ /T, where 
(RS) stands for Rees-Sciama. As a second step, we consider that there are Primordial fluctuations AT^ /T 
present on the LSS that will also be affected by the presence of the Void: in fact, they will be lensed by it, 
leading to an additional effect, AT^/T. So, finally, the total temperature fluctuation will be given by: 

AT Ar' p ' AT^ 3 ^ AT^' 

rj-i rj-i rj-i rj^ ' ( ) 



where by definition each (i)-labclcd fluctuation is given by 

T = f 



(2) 



with the bar representing the angular average over the sky and T = JV TW = 2.73A". 

At this point it is important to be more precise about the origin of the Void, relatively to the Primordial 
fluctuations. If such a structure comes from the same physical process as the Primordial Gaussian spectrum, 
then this would mean that there is some correlation between the two. However we will assume, from now on, 
that the location of the Void in the sky is not correlated at all with the Primordial temperature fluctuations 
coming from inflation. This is true, for example, if such structures come from a different process, such as 
nuclcation of bubbles. Note that this is a conservative assumption since, in the presence of a correlation, 
the temperature correlation functions would be generically enhanced. 

Then, we need to fix the properties of the Void in order to compute the RS (and Lensing) contributions. 
Physically, one may characterize such a Void via the following quantities: its comoving distance D from 
us, its comoving radius L, its present-day density contrast at the centre S$, and some shape for the density 
profile. As we will discuss, the RS fluctuation AT^^ /T turns out to be described by two parameters: its 
amplitude at the centre A and its angular extension, characterized by the diameter a of the Cold region. 
We also need some shape for the temperature profile, which will be determined by choosing some shape 
for the density profile of the Void. Let us discuss in more detail how the parameters A and a arc related 
to the physical parameters of the Void, namely D,L,Sq: the amplitude A mostly depends on the radius 
L and the density contrast 8q, while it is is not very sensitive to D; the diameter a depends on the ratio 
L/D. As explained in the next paragraph, we choose the values for A and a phenomenologically. Clearly 
this leaves a degeneracy in the choice of the physical parameters, since we have three of them (D, L and 
Sq) and only two observational constraints (A and a). 

In order to fix a range for the amplitude A we rely on the values given by|23[ as follows. The minimal 
temperature in real space observed in the Cold Spot area is AT ~ — 400/iK [J], but this value includes as 
well a contribution from the Primordial Gaussian profile, AT' P '/T. The authors of [H use a temperature 
profile for the secondary effect, add the Gaussian fluctuation and fit the observational data: this leads 
to an estimation for the value of the secondary effect at the centre, given by AT = —(190 ± 80)/iA, or 
equivalently A = (7 ± 3) x I0 -5 . In this paper we assume the latter range of values for A. This is probably 
not entirely accurate (the shape of the profile that we use is slightly different from theirs, since we consider 
a compensated Void, while the profile in [23] corresponds to the profile due to a cosmic texture), but it 
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should give a good estimate for the range of values of interest. In any case, our results can be easily rescaled 
for different values of A. As for the angular size a of the profile, we proceed as follows. According to [EHI, 
a shape that fits well the Cold Spot has a diameter of about 10° for the very cold part, but we also show 
in all plots the extreme case a = 18°: in fact as one can see in fig.l of [(| the entire cold region extends up 
to such large size. 

Having fixed the numerical values, we choose a metric to model such an inhomogencous region. The choice 
that we make is a spherically symmetric LTB metric, which is matched to a homogeneous and isotropic 
Friedman-Lemaitre-Robertson- Walker (FLRW) flat model. For simplicity we consider a flat FLRW Uni- 
verse with f2j\f = 1- Including a dark energy component would lead to an ISW effect already at the linear 
level, in addition to the RS effect that we consider, which is present in any case, even in absence of dark 
energy. The effect of a cosmological constant in an analogous setup has been studied e.g. by [1, [1] and the 
two effects turn out to be similar, with the same dependence on the Void radius, but with some difference 
on the dependence on the density contrast So 2 ■ 

Another important point is that our profile is a compensated Void, with the undcrdense central region 
surrounded by a thinner external overdense region. This is a requirement dictated by the matching condi- 
tions, which corresponds to the condition that the Void does not distort the outer FLRW metric. Physically, 
if the Void comes e.g. from a primordial bubble of true vacuum, this is a consistent requirement: in fact a 
bubble would have a thin wall with localized gradient energy, compensating for the lower energy contained 
in the true vacuum in the interior region. Note that, because of this feature, the angular size of the entire 
LTB patch that we consider will be larger than er, since we have a hot region as well. In any case, we 
show that the main contribution to the power-spectrum and bispectrum comes from the inner undcrdense 
region. 

Given these specifications for AT^^/T we now compute its shape in the following subsections. 



A. Temperature Profile 



First of all, we model the profile as an LTB metric with irrotational dust, which describes a compensated 
Void. It has been shown in [24| that this metric can be treated, in some cases, as a perturbation of an 
FLRW metric with a gravitational potential $ given by the following (Newtonian gauge) expression 

ds 2 = a 2 (r) [-(1 + 2$)dr 2 + (1 - 2<S>)dx l dx j ] , (3) 

where r is the conformal time and x 1 are dimensionless comoving coordinates. We have chosen units such 
that the present value of the conformal time is To = (67r) -1 / 6 . As a function of the dimensionless comoving 
radial coordinate r, the gravitational potential is given by: 

931/3 rr 

Hr) = —z vaTs / k{r)rdr , (4) 



5 (2tt) 2 / 3 

where k(r) is an arbitrary function which represents the local curvature and determines the shape of the 
density profile. This approximation is valid as long as k(r) is small. The only constraints on this function 
comes from the smoothness of the density profile at the centre (which dictates k'{0) = 0, where the prime 
denotes a derivative with respect to the r coordinate) and from the requirement that the LTB patch 
matches to a flat FLRW universe (fc'(L) = k{L) = 0). We have chosen arbitrarily the function k(r) as 
follows: 



k(r) = fc 







1 - 









^ = 2W76^°< (5) 



2 Specifically, the 80 dependence is different in the presence of A: the second-order term (RS effect) is suppressed, but there 
is a non-zero linear term. For small density contrast, 8q -C 1, the linear ISW would dominate in a A-dominatcd cosmology. 
The quadratic term dominates if So is large or if Qm is close to 1. An interesting case is the one for intermediate values 
(<5rj fa 0.3): the two effects are roughly compensating and the calculation does not depend much on the value of Qm- 
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FIG. 1: Plot of the present day density contrast 5(r), normalized for a central value of So = —0.2. The solid 
(dashed) line corresponds to a = 4 (a = 2). 



where a > 1, tl is the dimensionless radius of the LTB patch and Hq is the present value of the Hubble 
parameter. The density contrast that results from the choice of k(r) is approximately given by the following 
expression [HI]: 

where p(r) is the matter energy density inside the LTB patch, p is the average density in the outer FLRW 
region. Note that the density contrast can be large even if k(r) is small. We show in fig. [T]thc density 
profiles corresponding to certain values of a. From ([5]) and © it follows that the value of the normalization 
constant ko is directly related to So (the value of the present day density contrast at the centre of the Void): 

l _ 20ir5„ 
ft - 3(l-do) • 

Given the potential <£>(r) we can compute the AT/T for a photon which travels through this patch, 
by computing a line integral, following the expression given in [25l |2(| which is valid at second order in 
perturbation theory in $: 

Arp(RS) PT E /-, in \ 

f = 2 /,/"(^-^°)- <7 > 

where te is the value of the conformal time at emission (although of course the integrand is nonzero only in 
the region where the LTB structure is located). The conformal time evolves simply as r(r) = (r — ro) + To, 
where the subscript O refers to the observer space-time point. Here To is a non-local quantity, which is 
given in terms of <I> as [13, HE] : 

T = - T % f dr [$(r)' 2 +2f*'(f)*"(f)] . (8) 

Then, using the spherical coordinates angles, < 6 < tt and < <fi < 2n, we may parameterize the 
temperature profile as 

AT {RS \„ f A f{6) if6<e L „ L 

(e,6) = { Jy> L , tan6» L = — , (9) 

T v ' v> I if 6 > 6 L ' D ' w 



where the profile function f(0) is normalized so that /(0) = —1. The profile has no dependence on </> because 
we choose the z axis to point towards the centre of the Void. The A factor can be computed performing 
the integral ([7J analytically along a radial trajectory, and it is a function of the physical parameters of the 
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FIG. 2: Plot of the profile f(9) for 9 L = 20°. The dark solid (dashed) line corresponds to a = 4 (a = 2) and a 
Cold Spot with diameter a — 18° (<j = 16°). The light solid line shows, for comparison, the temperature profile 
(normalized so that AT/T — — 1 at the centre) of the texture that gives the best fit to the Cold Spot, as claimed 
in [13. 



configuration: the radius L of the LTB patch, the density contrast at the centre of the Void So and the 
distance D between the observer and the centre of the Void. The result is 3 , at leading order in LHq, 

Note that the dependence on the distance D is a weak correction, unless the patch happens to be located 
at distances comparable to the horizon. We can compute the function f(9) by numerical integration of ([7|) 
along a non-radial trajectory (computed as an unperturbed straight line, since any deviation would lead 
only to higher order corrections). In fig. [2] we plot the profile f(0) as a function of 9 for two values of a. 
For the reader's convenience, we provide polynomial interpolations of the profiles: 

f{6) = -1 + 6.663a; 2 - 5.954a; 4 - 17.258.x 6 + 33.959a; 8 - 19.361a; 10 + 2.940a; 12 , for a = 4 

(11) 

f(6) = -1 + 11.191a; 2 - 37.576x 4 + 58.272a; 6 -46.190.t 8 + 17.904a; 10 -2.601.T 12 , for a = 2 , 

where x = 6/6l- We have added also a light solid line in fig.[2]which shows, for comparison, the temperature 
profile of the cosmic texture that provides the best fit to the Cold Spot, as discussed in [23[ (its profile has 
also been normalized so that AT/T = —1 at the centre). Notice that the compensated Void has a hot ring 
in the profile, while the texture does not have it 4 . 



3 The dependence on the shape of the density profile (the a parameter) comes in the numerical factor and it is given by 

_/„ ■ — _„ o 10 , 4 " , — — r^v, which ranges between 0.05 and 1.85. This variability in the prefactor could account for the 

7(8q 4 + 50q 3 + 105q 2 +90c<+27) ' ° j r 

prefactors obtained by previous analyses @,I3i which are of about 0.1. 

4 There is one caveat, however: in a A dominated cosmology and with small contrast <5o, the linear ISW can become the 
dominant effect; and it has been shown by Q that this case does not lead to a hot ring. 
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FIG. 3: Plots of a l0 /A as a function of i. Left: The dark (light) solid line corresponds to a profile with a = 4 
and 9l — 20° (11°), namely a Cold Spot with diameter a = 18° (10°); dashed lines correspond to similar profiles 
but with a = 2 and a — 16° (9°). Right: The solid line corresponds to the full compensated profile with a — 4 and 
8l = 20°, while the dashed line to the non-compensated profile, namely the previous one truncated so to include 
only the cold region (Void). The light curve shows, for comparison, the ratio aeo/A that one would obtain from the 
temperature profile of the texture giving the best fit to the Cold Spot [53] . 



B. Decomposition into spherical harmonics 

The spherical harmonic decomposition of the (i)-th profile for the temperature anisotropy AT' ! '(n)/T 
of |T]) is defined as: 

41 ^Jdn r; m (n) . (12) 

For the (RS) component, since the profile is symmetric with respect to the z axis pointing towards the 
centre of the Void, only the ai m coefficients with m = are non- vanishing (and they are real). In the left 
plot of fig. [3] we show the ratio /A as a function of the multipole £. The dark solid (dashed) curve 
corresponds to a temperature profile with 9l = 20° and a = 4 (a = 2), hence to a Cold Spot with diameter 
a = 18° (cj = 16°); the light solid (dashed) curve is obtained by choosing instead 9l = 11°, so that a = 10° 
{a = 9°). Notice that the medium amplitude of such a ratio is roughly equal to the fraction of the sky 
covered by our LTB patch, namely about 3% for 6l = 20° and 1% for 6l = 11°, respectively concentrated 
at multipoles 10 < t < 20 and 15 < I < 40. 

One may wonder whether the presence of a compensating hot shell (see fig. [2]) in our profile has a 

significant impact on the magnitude and shape of the a eo coefficients. As an example, in the right plot 

of fig. [3] we focus on the with 0l = 20° and a = 4 (solid line), and show how they would change 

by truncating the temperature profile only to the cold part (dashed line). This approximately mimics a 
truncation of the density profile only to the underdensc region (the true Void) , disregarding the overdense 

compensating shell. One can sec that the a i0 are now suppressed but still quite similar, so that for our 
purposes the result is not going to be very different. This shows that the RS effect mainly comes from the 
underdense region. The same plot also shows, for comparison, the ratio aio/A that one would obtain from 
the temperature profile of the texture giving the best fit to the Cold Spot [23|. Notice that the texture 
and the non-compensated Void are pretty similar. 
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III. POWER SPECTRUM 



Given a temperature profile with its ag m coefficients one can compute the associated two-point correlation 
function. In general, given a single set of ag m coefficients, the two-point correlation function (power 
spectrum) is defined via the Cgs coefficients as 



c^ESttt- as) 



2^ + 1 

Note that this definition ensures that the C/s do not depend on the choice of the coordinate system on 
the sphere. Therefore we are free to keep our z axis aligned with the centre of the Void. In our case 
aim = a inf + 4m and, in order to face it with the experimentally observed value of Cg, we 

have to estimate the theoretical prediction for (Cg), where the (...) brackets stand for a statistical average 
over an ensemble of possible realizations of the Universe - or, equivalently, an average over many distant 
uncorrelated observers. 

For a Primordial and Gaussian signal the two-point correlation functions are given by: 

(^ 1 «a;)=^3^ im2 (C ) ). (1 4 ) 

(P) 

where the (C\ ) are predicted by some mechanism {i.e. inflation) that can generate Primordial Gaussian 

fluctuations. Then, there are two types of effects on the power spectrum due to our secondary effects: 

( p) 

i) the inflationary prediction Cg = {C\ ) receives corrections; 

ii) there is also some non-diagonal correlation between different £'s. 

As we have already stressed, we assume that the Rees-Sciama (RS) component in {TJ is uncorrelated 
with the Primordial (P) component, which means that it can be factored out of the brackets. Given this 
fact, from |T]) and (fl~2|) . we get two types of non-zero contributions for (Ce): the RS-RS contribution 

proportional to \a\ m | , and the P-L contribution proportional to (a) J a\ m ). The first effect is computed 
in the next subsection, while the second in the companion paper [14(. 



A. Rees-Sciama power spectrum 

As we have seen, for a spherical Void in the ^-direction, a,g m = if m ^ 0. Hence the RS contribution 
to the power spectrum coefficients, (Cg) = (C^) + Cj fi , is just: 

r (RS) _ \a l0 h ,s 

In the left and right panels of fig. 0] we have plotted respectively the quantities C^ RS ^ an( j 
{Cg) 1 v~ 1 ' Tq , for the range of parameters mentioned in sect. HU The shaded regions are indeed obtained 
varying the amplitude A in its l-cr range, namely A = (7 ± 3) x 1CP 5 . The results for different values of 

A can be extracted just by noticing that C^ R oc A 2 . As one can see, for a Void with a size which can 
account for the Cold Spot, the correction to the power spectrum is non-zero only in the range 5 < t < 50 
and its magnitude is about 5% — 25% with respect to the Primordial signal. We also note that the RS 
correction is of the order of the cosmic variance, ACf = (C ( e P) )^2/(2£+l). For comparison, we also show 
in the left panel of fig. [4] the correction to the power spectrum obtained from the texture of (23[, which 
turns out to be smaller and shifted at lower multipoles with respect to the Void case. 

In order to see how large is the impact on the cosmological parameter estimation due the RS effect from 
a large Void, one should perform a detailed statistical analysis of the CMB data, which we postpone for 
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FIG. 4: In the left panel we plot the C^ RS ^'s from the Rees-Sciama effect. The dark blue (light green) shaded 
region corresponds to an angular diameter for the Cold Spot equal to a = 18° (10°), for a range of amplitudes 
A = (7 ± 3) x 10~ J . For comparison, we also show (light orange) the analogous correction obtained in the case 
that the Cold Spot is due to a texture [2^. In the right panel the RS correction is added to the best-fit Primordial 
spectrum, from the usual ACDM concordance model; the cosmic variance (gray band) and the experimental binned 
data points are also shown. 



future work |27| . For the purpose of this paper, we may just estimate the impact that the inclusion of the 
RS effect has on the \ 2 °f the fit, which gets increased by the amount: 



C (RS) 2 



2<i 

where a\ is the variance of the two-point correlation function, including the cosmic variance and the 
instrumental noise. For an experiment like WMAP, the cosmic variance is the dominant source of error 
at low-^'s, so that we can just neglect the instrumental noise and use a\ = (cj ) 2 2jTI- The sum gives 
roughly the following result, for different angular diameters a of the Void: 

for a = 10° 

for a = 18° V ' 

Finally, as already mentioned, there is also a non-diagonal contribution to the two-point correlation 
function, which correlates different £'s in the range 5 — 50 and is of the same order of magnitude as the 
diagonal contribution to the power spectrum, Ci R . 




IV. BISPECTRUM 



Having computed the ae m coefficients and the two-point correlation functions, we now estimate the 
impact that a large Void has on the bispectrum coefficients, also to check whether this observable could 
be used to constrain the size and the density of the Void itself, which is a highly non-Gaussian object. 
Moreover, since the bispectrum is the main tool used to make detections of a primordial non-gaussianity, 
it is interesting to see whether the presence of a large Void can contaminate this detection and to what 
extent. 

Note that the bispectrum from the the RS effect due to conventional structures, i.e. with non-linearities 
at scales of order 10Mpc/ft., has already been studied by several authors [28|, 29]. However, the calculation 
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here is quite different for two reasons: i) we deal with much larger values for the a^ m 's, because we consider 
larger values of L and the RS effect on the temperature profile is proportional to L 3 ; ii) we are assuming no 
correlations between the RS and the Primordial fluctuations, while this is not the case in the conventional 
calculation where the RS effect arises from the structures seeded by the Primordial fluctuations themselves. 
In addition, we compute in [l4j the associated Lensing effect. 

The basic quantities are now the -B™ 1 ^ 2 ™ 3 coefficients, defined as 

ti t 3 is = "'lmi tt «2"i2 a i 3 m 3 , {i-O) 

which are coordinate-dependent quantities. So, in analogy with the C^'s coefficients, one introduces the 
angularly averaged bispectrum 



B *Ms - 1^ [ mi m , m , j B h i 2 i 3 > ( 19 ) 

mi,m2 ,772.3 



mi m.2 1713 



where the matrix represents the Wigner 3-j symbols and the sum is carried over all possible values for the 
m^'s. One can indeed show, by using the properties of the Wigner 3-j symbols [30j |. that this combination 
does not depend on the chosen z-axis, so that these quantities are more suitable to make predictions. For 
convenience, we nevertheless keep our z-axis along the direction of the centre of the Void. 

We are thus interested in evaluating (B^ 2 ^ 3 ) and, using ([I]) and (|19j) . one realizes that it corresponds to 
a sum of 27 terms, of which 23 have zero statistical average. As already mentioned, a crucial assumption 
that we make here is that the coefficients not stochastic quantities, which means that the location 

of the Void in the sky is not correlated at all with the Primordial temperature fluctuations coming from 
inflation (this is a conservative assumption: if there is some correlation, the non-gaussianity could be much 
more important, since some terms would be non-zero). Under this assumption, the four types of terms that 
potentially survive are the ones involving: ((a ( - RS) ) 3 ), (a^a^a^ 5 )), (a ( - RS '> (a (i) ) 2 ) and (a^^a^) 2 }. 
However, as shown in Appendix |A"1 the potentially very large terms (a^ 5 -* (a^ p ^) 2 ) are actually exactly zero 
(because of the absence of correlations between RS and P). Moreover, since a^ L ' <C a^ p \ we can neglect 
the (a,( RS ' i (a( L )) 2 } terms with respect to the (a^^a^'a^ 5 '} terms. 

Summarizing, we are left with two types of contributions to {B^g^): 

/R (i?S) \ \- ( t\ h 4 A /(RS)JRS) (RS), (t) „, 

\ D ixt?,t 3 l - 2^ I mi m 2 TO3 / * imi <2m 2 %m 3 / ' \ zu > 

mi,m 2l m 3 

We compute the first one in the following subsection, while the second in [T3 |. 



A. Non-gaussianity from RS effect 

To compute the (Bgff 2 ] 3 ) term is very simple. The only non-zero a^^f^ coefficients are those with m = 0, 
so that: 

b (rs) _ ( ii £2 4 A as RS rs (22) 
n eM 3 ~ \ / tl0 i2 ° ^° ' ^ ' 

It is customary to define a reduced bispectrum bg^g^ via the following: 



(24 + l)(2f 2 + l)(24 + l) li h I 



Btihh - \l ^ ( o o o ) b ei e 2 e 3 ■ (23) 
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FIG. 5: Plot of 10 16 ^V + 1) 2 /(2tt) 2 b ( e f e S) as a function of the multipole i, for a = 18° and a = 10°, with 
A = (7 ± 3) x 10" 5 . For comparison, we also plot in red the prediction for the primordial signal for /ml = 10 . 
The light (orange) shadow shows, for comparison, the result for the texture considered in [23l ). 



We plot the diagonal contribution b iU in fig. [5J the equivalent amplitude is very high compared to a 
typical primordial signal, since it corresponds roughly to /jvz, ~ 10 4 . However, for i > 60 the RS signal due 
to the large Void goes rapidly to zero. Therefore, for experiments like WMAP or Planck, which go up to 
about I ~ 800 and I ~ 2000 respectively, only a small subset of the data is affected by the RS contribution. 

Focusing on the RS signal, we now turn to estimate whether it is detectable or not. The signal on a 
single multipole is lower than the cosmic variance: so we have to sum over all the Ps to find a bispectrum 
Signal-to-Noise (S/N) ratio. For a signal labeled by i, this is defined as (see e.g. [30|): 

F U = 2 ■ ( 24 ) 

where <Ji 1 t 2 t 3 is the variance of the bispectrum: 

Otitsta ~ (^2X^3)^1^3 ) ( 25 ) 

and A^ 1 £ 2 £ 3 = 1, 2, or 6 respectively if all ^'s arc different, if only two of them are equal or if they are 
all equal. The Ce's represent the sum of the CMB power spectrum plus the power spectrum of the noise 
of the detector. In general, at some i max the noise becomes dominant, while below i max the variance is 
dominated by the Primordial one, namely Ce — (C^). 

Neglecting the instrumental noise, in the left panel of fig. [H we show the result for (S/N) us as a function 
of imax- As one can see, the signal is detectable for a large part of the parameter space: so it should 
already be possible to look for such a signal in the WMAP data. Conversely, the absence of any detection 
would give interesting constraints on the size L and density 60 of a large Void. We show in the right 
panel of fig. E] the region of physical parameter space which would give rise to a detectable signal, namely 
(S/N)rs > 1: for a cold region diameter a > 10°, this happens for all the points in the light shaded 
region. For a > 18° such region gets larger and includes the dark shaded one. For instance, for a Void 
with cr = 18° and So = —50%, a signal in the bispectrum would appear only if L > 200Mpc//i, which 
corresponds to A > 4 x 10~ 5 . 

Other candidates for a structure that could explain the Cold Spot might be subject to constraints 
analogous to those discussed above for a big Void. For instance, from the left panel of fig. [5] one can see 
that the texture proposed in [23[ would give rise to S/N > 1 in the bispectrum only if A > 7 x 10~ 5 . 



(S/N) t 



FZ 1 
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FIG. 6: Left: Plot of S/N for the RS as a function of the multipole £ m ax for a — 18° and a = 10°, with 
A = (7 ± 3) x 10~ 5 . The light (orange) shadow shows, for comparison, the analogous quantity for the texture 
considered in [2^|. Right: The black lines, from bottom to top, show the contour levels for A — (4,7, 10) x 10~ 5 
in the plane (L, \So\). For a = 10°, all the points of the light shaded (green) region would give (S/N)rs > 1- For 
a = 18° this region gets larger and includes the dark shaded one, which goes down to A = 4 x 10~ 5 . 



B. Contamination of Jnl measurements 



We now turn to the impact that a huge Void in the line of sight would have on the measurement 
of the primordial non-gaussianity parameter /nl, in terms of which it is customary to parametrize a 
primordial non-Gaussian signal. By definition, /jvl is introduced (see e.g. [3(j for details) parameterizing 
the primordial curvature perturbations 4>{x) as follows: 

4>{x) = 4> L (x) + fNL^l(x) - (^(x)» (26) 

where 4>l(x) is the linear Gaussian part of the perturbation. Given a physical model (like slow-roll inflation) 
/jvx is generically a function of the momenta, i.e. fNL^k), but in the quantitative data analyses it is usually 
assumed to be just a constant number. Note that single field minimally coupled slow-roll inflationary models 
predict very small values for Jnl, that is /jvl = O(0. 1) [Ml, [11], but other models may predict larger values 
(see e.g. [301]). The primordial bispectrum coefficients can be written as: 

BfZl=fN L Bf^ 3 , (27) 

where the Bf™™ have a specific form in terms of the primordial spectrum and the radiation transfer 
function. 

As we have seen, the RS effect leads to a large contribution to {Bg 1 g 2 i i ) for multipolcs in the range 
10 < i < 50. On the contrary, the Lensing effect [I3| is much smaller at low £'s, but could contaminate 
the primordial bispectrum signal at large £'s, since it couples the low RS-^'s with the high ts of the 
primordial signal 11411 . The impact on /jvl due to the RS effect can be computed by estimating the 
following ratio [30l. |33| : 

Af N T(£ max )= '±\ ■ (28) 



E 



2<t t <l 2 <t 3 <l maI 



In other words, if a large Void exists, one should subtract it from the data in order to get the correct value 
for /jvl, thus avoiding to overestimate the latter by the amount Af^f . 
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FIG. 7: The correction to /jvz, arising from the inclusion of the RS effect, A/^^ (£ m ax), plotted as a function 
of £max, for a = 18° and a = 10°, together with A = (7 ± 3) x 10~ 5 . Left: the log-log scale allows a better 
understanding of the effect for different experimental sensitivities. Right: the linear scale is more useful for a 
comparison with the constraints on /jvz, reported by the WMAP 1-year analysis [35| . The light (orange) region 
shows the analogous results in the texture case [231 ]. 



We may very easily give an approximation of the 2?^™^ appearing in the numerator of (|28[) by using 
with the Sachs- Wolfe approximation for the primordial signal: 

c^a = -6«cr>< c r> + (^fxo + (^rxO) • w 

This is a good approximation only for the low-f plateau in the power spectrum, while the full expression 
should be used for higher ts. However, since (-B^f] ) vanishes if any of the fs is larger than about 60, we 
consider it to be a sufficiently fair approximation to get an estimate of the numerator of l[2"8]) . This fact also 
allows to neglect the experimental noise in the numerator's <^f 1 g 2 i 3 - Instead the denominator is a quantity 
very sensitive to the experiment (see [30l]). In fact, at some high t max (dependent on the experiment) 
the experimental noise in the p\i 2 i 3 of the denominator becomes so large that the multipolcs £ > l max 
give a negligible contribution to the sum. Accordingly, the denominator of (|28p turns out to be equal to 
5.8 x 1CT 2 for WMAP (l max « 800) and 0.19 for Planck (l max « 2000) [H, so that one obtains: 



A/jf/M 1 fOTWMAP • (30) 
JNL '0.1 for Planck 



It is also interesting to study the dependency of &-f NL {(-max) for smaller values of l ma x- This can be 



easily done because the denominator in (|28p is well approximated by £ max x 10 [341 ] . until the experimental 
noise is negligible. The result is plotted in fig. [7] from the left panel it turns out that the RS effect due to 
the large Void does not affect much the /jvl measurements for high resolution experiments. In fact, the 
corrections are localized at low £'s because (Be 1 e 2 e 3 ) get a sizable RS correction only if £i, £ 2 , ^3 iS 60, while 
the search for a primordial non-gaussianity uses all the experimental data-points up to £ max ~ 800 — 2000. 
Therefore, even if the multipolcs £ < 60 are strongly affected, they only represent a small fraction of all the 
bispectrum data-points. Note also that the impact on /jvx is roughly the same for a = 18° and a = 10°: 
the reason is indeed that a smaller Void affects higher £'s, which are more relevant for the extraction of 
the primordial Jnl- 

The right plot in fig. [JJ allows a direct comparison with the WMAP-lyear experimental constraints on 
Inl [35[. We can already see that some region of the parameter space can be excluded for the Void with 
a = 18°: values for the amplitude A > 7(8) x 10~ 5 are excluded at l(2)-cr by the error bar localized 
at £ = 30. Of course, such an analysis would require a more refined technical treatment - for example 
including the recent WMAP 5-year data, sky-cuts, the full expression for B^J^ - which we do not address 
in the present paper. A full analysis should also include the correlation matrix with other cosmological 
sources: SZ-lcnsing effect, point sources and the primordial signal. Anyway, looking for /jvl is not the best 
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way to constrain a large Void: one should rather compare directly the observed bispectrum data with the 
prediction for (Bg R f\). Since we have shown in fig. E] that the RS Signal-to-Noise is above unity for most 
of the parameter space, a full analysis should be able to set interesting constraints. 

Finally, we may also give a comment on the results by (36j , whose authors claim a detection of primordial 
non-gaussianity already in the WMAP data. In their analysis they also find that, masking the region which 
corresponds to the Cold Spot, their measured /jvl goes up by an amount of 7. In the light of our analysis, 
this cannot be due to the RS effect, because the RS contribution to /jyx is smaller (about 1) for WMAP 
and, moreover, has opposite sign. The possibility that this could be ascribed to the Lensing effect will be 
discussed in 1141 1 . 



V. OTHER VOIDS 



We briefly comment in this section about the possibility that several Void regions are present in the sky. 
The authors of [HI, [TtJ have claimed detection of the ISW effect because it could explain the correlations 
that they were able to identify between the galaxy surveys and the CMB data. By exploring a region of the 
northern hemisphere that covers about 20% of the sky, has catalogued about 50 of such Voids 5 , which 
have a mean radius of about 5° and a mean amplitude A m 3.6 x 10 -6 (or, equivalently, AT sa — 11/iK); 
their density contrast for the luminous matter can be directly observed in the galaxy surveys; the knowledge 
of their redshift also allows to estimate their physical size L, which turns out to be of about 50 — 100Mpc//i. 
As stressed in the introduction, [l8[ has interestingly pointed out that the existence of such large Voids is 
unlikely to be explained by the standard structure formation scenario. It could be relevant, therefore, to 
apply our considerations to these Voids as well. 

If N Voids are present in the sky, however, there are differences with respect to the analysis developed in 
the previous sections. If we focus on the two-point correlation function, in addition to the N terms of the 
same type as (fT5"|) . we have to compute also N 2 — N interference terms. As for the three-point correlation 
function, it would contain N terms of the same kind as (f2"2"j) . plus N 3 — N interference terms. For a random 
distribution of Voids in the sky, one expects the interference terms to add up incoherently, leading only to 
an oscillatory modulation of the correlation functions. 

For a precise computation of the effect of these Voids on the power spectrum and bispectrum, it would 
be necessary to know, for each Void, its location, the amplitude A and the angular size a of its temperature 
profile. Note that the most of the contribution is expected to come from the few Voids with the largest A 
and a. The authors of [ijj give these informations for the 50 Voids identified, except for the amplitude A, 
for which only the mean value A is provided. 

We can nevertheless give an estimate of the impact on the CMB due to the RS effect of the Voids 
catalogued in [l7| (their location is sketched in the right panel of fig. [5]) , by assigning to each of them an 
amplitude equal to the mean one. The angular size of the Voids spans from 3° to 14°, and there are 19 
Voids with a > 10°. The left panel of fig. [8] shows the result for the two-point correlation function (where 
we also show that the interference terms are negligible) , whose order of magnitude can be understood by 
the following easy argument. Rescaling to the case of an amplitude A the results obtained in sect. IIIII for 
one Void with radius of about 5°, one gets: c\ RS ^ Tq ~ 0.1/xK 2 . Assuming that the interference terms 
(which depend on the relative location in the sky of the Voids) are irrelevant, the latter value for a single 

Void has to be multiplied by N = 50, thus obtaining C^' ^ +1 ^ Tq = C(5^K 2 ), consistently with the left 
plot of fig. [HI We may estimate the contribution to the three point function similarly. The signal-to-noise 
ratio for a Void of radius 5° and amplitude A would be of about 10 -3 , as shown in sect. HVI As done before, 
the final result is obtained by multiplying by the number of Voids, leading to (S/N)rs = 0(5 x 10~ 2 ), 
assuming again that the interference terms are negligible. 



5 The Cold Spot, located in the southern hemisphere, is clearly not among them. 
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FIG. 8: In the left panel we plot (solid line) the RS power spectrum coefficients obtained from the 50 Voids identified 
in pjl, for all of which we assigned the mean amplitude A = 3.6 x 1CP 6 . These Voids have a in the range 3° — 14°: 
the dashed curve is obtained by selecting (among the 50 Voids) the 19 ones having a > 10° . The dotted thin curve 
is the full result without interference terms. Right: a sketch of the location of the 50 Voids [l7| . 



VI. CONCLUSIONS 



Motivated by the so-called Cold Spot in the WMAP data, we have studied in this paper the impact on 
statistical CMB predictions, in particular the two and three point correlation functions, of the presence of 
an anomalously large Void along the line of sight. Indeed, the existence of such a Void could be at the 
origin of the Cold Spot, due to the Rees-Sciama effect (the Lensing effect is analyzed in the companion 
paper [H||). 

First, we have computed the temperature profile using an LTB solution of the Einstein equations, 
matched to an FLRW metric. Then, we have computed its impact on statistical predictions for the CMB, 
assuming this structure to be uncorrelated with the Primordial fluctuations. As suggested by [23j], we 
consider the angular size of such a Void to be about 10° — 18° and its temperature at the centre to be 
characterized by AT = -(190 ± 80)fj,K. 

For the power spectrum the results are as follows. The RS effect is non-negligible: we predict a bump 
of 5% — 25% to be added to the Primordial spectrum, localized at multipoles 5 < £ < 50. This should lead 
to a variation in the x 2 f° r the WMAP fits of order unity. 

Then we have studied the impact on the bispectrum coefficients. For the RS effect we have found that 
the Signal-to-Noise ratio is larger than unity at I > 40 for most of the parameter space, and therefore 
this should already be visible in the available data. Through the bispectrum, we have studied also the 
impact of such a structure on the determination of the primordial non-gaussianity. The RS bispectrum 
signal is large but localized at low multipoles (10 < I < 50), so it has a small impact on high resolution 
experiments, which can go up to very large multipoles: the overcstimation of f^L due to the RS effect turns 
out to be Af^f ) ~ 1 for WMAP and A/^ 5) ~ 0.1 for Planck. Using the already existing WMAPl-year 
constraints [351 ] on f^L at low £, one can exclude extreme values of the temperature contrast of the Void. 
For example, values for AT/T larger than 8 x 10~ 5 for a Cold Spot with diameter of 18° are likely to be 
excluded, via a full analysis. So, we conclude that the bispectrum is a valuable tool for constraining an 
anomalously large Void, through the RS effect. 

We have also considered the 50 Voids whose detection has been claimed in [ijj ■ These Voids have mean 
angular diameter of about 10° and average temperature at the centre characterized by AT w —ll[iK. The 
effect on the two-point correlation function is about 0.5%, hence much smaller than the one due to the 
large Void considered to explain the Cold Spot. The effect on the three-point correlation function has not 
been studied in detail but it is expected to lead to a Signal-to-Noise ratio smaller than one. 
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Finally, we have also applied our considerations to the case in which the Cold Spot is explained by a 
cosmic texture [23| rather than a large Void: the effect on the power spectrum is similar but somewhat 
smaller; the three-point correlation function leads also to a smaller Signal-to-Noisc ratio, which however 
turns out to be above unity for some part of the parameter space. 
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APPENDIX A: VANISHING OF THE P-P-RS CONTRIBUTION TO THE BISPECTRUM 

We show here that the term containing } in the bispectrum coefficients (B^i^) is van- 

ishing. More generally, any term of the kind ) is zero, for any a w uncorrected with a (p) . 

Averaging (JTHJ) over a statistical ensemble of possible realisation for the Universe one gets: 

{B m ilXT 3 ) iPPRS) = SmAaffo <m 2 af ams ) + (1, 2,3^2, 3, 1) + (1, 2,3^3, 1, 2) (Al) 
= <W a^R™ 3 « m2 <,- m3 > + (1,2,3^2,3,1) + (1, 2,3^3, 1, 2) 
= <W afJ(-r 3 C £2 ^ 2 , 3 5 m2 ,_ m3 + (l,2,3^2,3,l) + (l,2,3^3,l,2) . 

Summing over the m,'s according to (|19p . one finds 

(^ 1 Wa> (PPHS) =«f 1 o(^ 2 )^ A E(-) m2 ( m 2 _t 2 )+(l,2,3^2,3,l) + (l,2,3->3,l,2) = , 



(A2) 

where we have used a known property of the Wigner 3 — j symbol and the fact that by definition a p s = 0, 
see ©. 



